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<N ■ Abstract 

In the light of the recent LHC data on proton-proton and lead-lead collisions we examine the 
question of the multiplicity scaling of HBT radii in relativistic nuclei and particle interactions. 
Within the UrQMD transport approach we study a large variety of system sizes at different beam 
energies and extract the HBT radii. In the calculation, we find a good scaling of the radii as a 



function of charged particle multiplicity, if the change in the multiplicity is caused by a change of 



centrality at the same energy. However, the scaling is only approximate when the energy, y/s, is 
changed and breaks down when comparing pp to A A reactions. 



* graef@th.physik.uni-frankfurt.de 



I. INTRODUCTION 



The properties of strongly interacting matter are described by the theory of Quantum- 
Chromo-Dynamics (QCD). To explore the details of QCD matter under extreme conditions, 
one needs to compress and heat up QCD matter to regimes present microseconds after the 
Big Bang. Today these conditions can only be found in the interior of neutron stars or 
created in heavy-ion collisions at relativistic energies. Over the last decade the experimental 
programs at the SPS (e.g. with the NA49, CERES and NA50/NA60 experiments) and at 
RHIC (e.g. PHENIX, STAR, PHOBOS, and BRAHMS) have provided exciting pioneering 
data on the equation of state, the transport properties of the matter created and its spatial 
distributions |ll-llO|. These programs are currently extended into a system size scan with 
NA61 at SPS and a systematic beam energy scan with the RHIC-BES initiative. In addition, 
at the high energy frontier unprecedented data from the Large Hadron Collider (LHC) for 
(high multiplicity) proton-proton and Pb+Pb reactions up to ^snn = 7 TeV has become 



available (see 



11 



12j for HBT related results). Particle correlations, i.e. Hanbury-Brown 



Twiss correlation (HBT) or femtoscopy allow to gain deeper insights into the emission pat 



terns and coherence regions of the matter created 
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161 ]. One generally assumes that the 



observed HBT radii scale with the charged particle density (or number of participants) as 
the charged particle density should be a good proxy for the final state volume [17| • However, 
the interferometry volume may not only depend on multiplicity, but also on the initial size 



of the colliding system [18J . Indeed, one of the surprising LHC results concerns the scaling 
violation observed in pp reactions as compared to AA reactions at lower energies at the same 
charged particle density. In this paper, we want to explore the spatial structure of the source 
created in collisions of various heavy ions at different energies and centralities to shed light 
on the observed scaling violation when going from proton-proton to AA collisions at the 
LHC. Other investigations on the charged particle yield scaling can be found in p, I20l-l22|. 



Results for PbPb and pp reactions at the LHC within the same model can be cound in 
23,124 1 . 



II. MODEL AND HBT CALCULATION 



For the present study we employ the UrQMD 25|, |26| transport model in version 3.3 (for 



details of version 3.3 see 



27 



281) 



HBT results from UrQMD see 



n he model can be downloaded from |29J. For earlier 
30H33J. UrQMD is a microscopic non-equilibrium transport 
model. It models the space-time evolution of nucleus-nucleus collisions from the beginning 
of the collision until the kinetic freeze-out. Particles are produced via hard collisions, string 



excitation and fragmentation and via resonance excitation and decay. 



For the calculation of the HBT radii we use the pion freeze-out distribution from UrQMD. 
Then we calculate the HBT correlation function by [l6|, [l7| 



C(q, K) = 1 + / d 4 x cos(q ■ x) d(x, K) 



(1) 



where C is the correlation function, q is the four-momentum distance of the correlated 
particles, K = (pi + p2)/2 is the pair momentum, x is the particle separation four- vector 
and d is the normalized pion freeze-out separation distribution, which is an even function 
of x. For the analysis in this paper all values are taken in the pair longitudinal comoving 
system (LCMS). Since UrQMD generates a discrete set of freeze-out points, the integral in 
Eq. [TJis substituted by a sum. 

The HBT radii Rij are obtained by fitting the function 



C(q, K) = 1 + \(K)exp 



- £ <M;4(K) 



i,j=o,s,l 



(2) 



to the calculated three-dimensional correlation functions. For the analysis in this paper the 
correlation functions are fitted over a range |q;| < 800 MeV/c for proton-proton collisions, 
|q;| < 300 MeV/c for carbon-carbon collisions and |q;| < 150 MeV/c for all other collisions. 
The difference in the momentum ranges is motivated by the fact that the width of the peak 
in the correlation function gets broader for smaller systems. Thus, the fit range is bigger for 
proton-proton and carbon-carbon, than it is for lead-lead collisions. 
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FIG. 1. (Color online) The three HBT radii i? D Mt) -Rsi<fe and -R/ong as a function of the charged 
particle multiplicity at midrapidity, (dN c h/dr)) l ' s and fixed kx = 300 — 400 MeV. The lines with 
symbols are the simulation results. The gray triangles, the black circles, the red squares and the 
green crosses are for lead-lead collisions at yfs = 2760, 200, 130, 62.4 GeV (in the same order) at 
0-5%, 5-20%, 20-50% and 50-80% centrality for the different points. The pink crosses are results 
for carbon-carbon at \fs = 200 GeV for the same centrality classes and the beige diamonds show 

n 

results for various multiplicity classes from proton-proton collisions [23]. Blue circles and brown 
squares depict results for central copper-copper events at y/s = 200 GeV and central lead-lead 
collisions at Ei a b = 158 GeV. The green stars are experimental results for central gold and lead 



collisions at kx = 300 GeV/c taken from 1 



12] 



III. SCALING OF THE HBT RADII 



Fig. [T] shows the three HBT radii Rout, Rside and Ri ong as a function of the charged 
particle multiplicity at midrapidity (\r]\ < 1.2 for pp and \rj\ < 0.8 for all other classes), 
(dN c t l /dr]) lj ' 3 and fixed kx = 300 — 400 MeV. The lines with symbols are simulation results 
for lead-lead collisions at ^/s = 2760, 200, 130, 62.4 GeV for 0-5%, 5-20%, 20-50% and 



50-80% centrality, for carbon-carbon at \/s = 200 GeV in the same centrality classes, for 
proton-proton at y/s = 7 TeV with different dN c h/drj classes, for central copper-copper 



ead collisions at Ei a b = 158 GeV. The 
For nucleus-nucleus reactions one 



ij. 



collisions at \fs = 200 GeV and for central lead- 
green stars are experimental results taken from 1 
observes a rather linear scaling with (dNch/dr]) 1 ^ 3 . The scaling is very good if the change in 
(dN ch /driy/ 3 is caused by a change of centrality at a fixed energy. However, a small offset 
on the order of 2 fm - 3 fm is visible for different system sizes, if the radii are extrapolated 
to N c h — > 0. This is expected due to the finite size of the nuclei in AA reactions [18}]. In 
contrast, increasing the center-of-mass energy leads to a reduction of the radii at a given 
fixed iV c /j-bin. The scaling of the source size with (dNch/dr]) 1 / 3 for different centralities is 
a hint that the underlying physics, e.g. pion production via resonance decay versus pro- 
duction via string fragmentation, is nearly unchanged by changes in the collision geometry. 
A change in ~^fs on the other hand results not only in different weights of the production 
mechanisms, but also in changed expansion dynamics towards a more violent expansion 
with increased energy. Qualitatively, one expects a scaling of the length of homogeneity 



as R = Rgeom/ y 1 + (vj_)7n±/2T 18|, Il9| . where R ge om is the geometric size of the collision 
region, v± is the transverse flow velocity and T is the freeze-out temperature. I.e. the 
increase in transverse flow leads to a decrease of the observed radii with increasing energy 
as observed in the model. This combination leads to a deviation from the (dN c h/dri) 1 / 3 
scaling of the HBT radii. The proton-proton calculation (and the data) show significantly 
smaller radii and a different slope from what is expected from nucleus-nucleus results. This 
behaviour is attributed to the strongly different particle production mechanisms in AA and 
pp. I.e., bulk emission vs. string/jet dominated emission which is also in line with the 
theoretically observed dependence of the HBT radii on the formation time of the hadrons 



from the jet fragmentation and string decay 23]. 



Since the K± dependence of the HBT radii tells us much about the expansion of the source 
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, ll5j], let us next investigate how a variation of dN c h/dr] is reflected in the differential 



HBT radii as recently discussed in |34{. Fig. [2] shows the three HBT radii R ou t, R S ide and 
Riong at fixed charged particle multiplicity at midrapidity as a function of kx- The shown 
calculations are chosen so that they fall roughly into two {dN^/dr]) classes. The first class 
contains calculations with {dN ch /drf) ~ 600 (exact values are 670 for Pb+Pb at y/s = 2760 



GeV, 20-50% centrality and 665, 595 and 509 for Pb+Pb at y/s = 200, 130, 62.4 GeV, 0-5% 
centrality). The second class contains calculations for (dN c h/drj) « 25 (exact values are 23 
for C+C at v/i = 200 GeV, 0-5% centrality and 32, 28 and 23 for Pb+Pb at y/s = 200, 130, 
62.4 GeV and 50-80% centrality). 

A very similar slope in K^ is observed for all UrQMD results. This leads to the conclusion 
that the observed HBT radii dependence on the radial flow in the model is weaker than 
observed in the data. The shift in magnitude of the radii is related to the magnitude 
differences already observed in Fig. [TJ that are mainly dominated by geometry and y/s 
effects. 

IV. VOLUME AND FREEZE-OUT TIME 

Next, let us investigate the energy and system size dependence of the homogeneity vol- 
ume. Fig. [3] shows the volume of homogeneity as a function of dN c h/dr) for various systems. 
Lead-lead calculations are shown for y/s = 2760, 200, 130, 62.4 GeV (grey triangles, black 
circles, red squares, green crosses) in the centrality classes 0-5%, 5-20%, 20-40% and 40-80%. 
The pink crosses show y/s = 200 GeV carbon-carbon results for the same centralities, and the 
beige diamonds represent proton-proton calculations at y/s = 7 TeV for different dN c h/drj 
bins. Blue circles and brown squares depict results for central copper-copper events at y/s = 
200 GeV and central lead-lead events at E Lab = 158 GeV. These results are compared to 
experimental data [ll-ll2| which is represented by green stars. In line with the experimental 
data, a strong increase in the volume proportional to the charged particle multiplicity is 
observed. A good agreement between experiment and theory is observed for the quantity 
Rl ide Ri ong while the experimental results for R out R S ideRion g are slightly overestimated. This 
is due to a too large R out in the calculations. The overestimation of R out is common for 



hadronic cascade mode 
heavy ion collision 31 



s and can be explained by a lack of pressure in the early stage of the 



35| . While the volume of the homogeneity region for each individual 
energy scales very well with dN c h/drj Fig. [3] shows a steeper slope with decreasing energy. 
The calculations also hint to an offset for AA reactions on the order of 25 fm 3 (i?^ ide i?/ onp ) 
and 50 fm 3 (R out R side R long ) . 

Finally, we explore the apparent freeze-out times r/. The results are obtained by fitting 
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FIG. 2. (Color online) The k± dependence of R ou t, Rside and Riong- The black dots are calculations 
at \fs = 2760 GeV and 20-50% centrality, the red squares, the green crosses and the pink crosses 
are lead-lead for 0-5% centrality at </s = 200, 130, 62.4 GeV. They have (dN ch /dr]} w 670, 665, 
595, 509. The other presented calculations are carbon-carbon at yfs = 200 GeV for 0-5% centrality 
(blue circles) and lead-lead at yfs = 200, 130, 62.4 GeV (beige diamonds, grey triangles, blue 
triangles) all for 50-80% centrality. These collisions have {dN^/drj) ~ 23, 32, 28, 23. The green 
stars represent ALICE lead-lead data for central collisions at y/s = 2760 GeV [12] . The blue 
diamonds are experimental results for central gold-gold collisions at y/s = 200 GeV from the 
STAR collaborator. |lj 



the hydrodynamically motivated Eq. [3] [12|, |36| to the k± dependence of Ri ong in the interval 
Kj_ = 200-800 MeV/c. For this purpose the pion freeze-out temperature is assumed to be 

T = 120 MeV. 

2 T K 2 (m ± /T) 



n 2 

■^long 



7 



m_|_ Ki(m±/T) 



(3) 



where m± = Jm 2 + k\ and Ki are the integer order modified Bessel functions. Fig. H] shows 
the freeze-out time as a function of dN ch /drj for various systems. The grey triangles, the 
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FIG. 3. (Color online) Two definitions of the volume of homogeneity as a function of energy for 
various systems. In the left plot the volume is defined as RoutRsideRiong an d in the right plot the 
volume is defined as R 2 side Ri ong . The gray triangles, black circles, red squares and green crosses 
depict UrQMD results for lead-lead collisions at (in this order) yfs =2760, 200, 130, 62.4 GeV for 
the centralities 0-5%, 5-20%, 20-40%, 40-80%. The pink crosses are carbon-carbon calculations at 
yfs = 200 GeV for the same centralities, the blue circles are central copper-copper collisions at 
y/s = 200 GeV and the brown squares are central lead-lead collisions at Ei a b = 158 AGeV. The 
beige diamonds depict proton-proton results at y/s = 7 TeV for different (dN^/dr]) 1 ' 3 classes. The 



green stars show experimental results taken from 



u 



12|- 



black circles, the red squares and the green crosses are calculations of lead-lead collisions 
at v/i = 2760, 200, 130, 62.4 GeV (in the same order) for the centralities 0-5%, 5-20%, 
20-40%, 40-80%. The pink crosses are carbon-carbon collisions at \fs = 200 GeV for the 
same centralities. The blue circles are calculations for central copper-copper collisions at 
Ws = 200 GeV and central lead-lead collisions at E^ = 158 AGeV. Experimental results 
[lldlOl Il2j are depicted by green stars. As for all the other observables, there is scaling for 
each energy individually. As anticipated from the calculations of R\ ong the decoupling time 
Tf increases with decreasing energy. This confirms the idea of a shorter decoupling time 
with increased energy. The offset in Tf for dN c h/dr) — > seems to hint towards a minimal 
decoupling time T} nm ~ 4 — 8fm/c in AA reactions and Tl nm < 2 fm/c in pp. 
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FIG. 4. (Color online) The freeze-out time as a function of energy for various systems. The gray 
triangles, black circles, red squares and green crosses depict UrQMD results for lead-lead collisions 
at (in this order) ^ =2760, 200, 130 62.4 GeV for the centralities 0-5%, 5-20%, 20-40%, 40-80%. 
The pink crosses are carbon-carbon calculations at \fs = 200 GeV for the same centralities, the 
blue circles are central copper-copper collisions at \fs = 200 GeV and the brown squares are central 
lead- lead collisions at Ei a b = 158 AGeV. The beige diamonds depict proton-proton results at y^s = 
7 TeV for different (dN^/dri) 1 ' 3 classes. The green stars show experimental results taken from 

OUQ. 



V. SUMMARY AND OUTLOOK 



In the light of recent LHC data on pp and AA collisions, which indicate a modification 
of the multiplicity scaling of the HBT radii, we have explored the A^ scaling for a large 
variety of systems and energies. We find good scaling of the radii with dN c h/dr) within a 
given system and energy. While the radii decrease slightly with increasing beam energy, 
they have a similar slope when plotted versus (dN^/dr/) 1 / 3 at all energies. When analyzing 
the freeze-out volume versus dN c h/dr] the increasing steepness of the slope for decreasing 
energies becomes visible. For all observables the scaling of the results for pp collisions differ 
strongly from the nucleus-nucleus results. We relate this observation to the different particle 
emission patterns (bulk vs. strings) in A A and pp. 
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